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We present results for the renormalized quartic self-coupling An and the renormalized Yukawa coupling i/r in 
a fermion-Higgs model with two SU(2) doublets, indicating that these couplings are not very strong. 
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1. INTRODUCTION 

It is an important issue to investigate within a 
non-perturbative regularization scheme whether 
the quartic self-coupling and Yukawa coupling of 
the fermion-Higgs sector of the Standard model 
(StM) remain relatively small when increasing the 
bare couplings to very large values. For this it 
is desirable to construct a lattice fermion-Higgs 
model with a realistic fermion content. A naive 
transcription of the continuum lagrangian with 
one SU(2) doublet leads to the large number of 
16 doublets on the lattice because of the species 
doubling phenomenon. There are two propos- 
als which allow to reduce this large number of 
mass-degenerate SU(2) doublets to one: The mir- 
ror fermion model 0| and the reduced staggered 
fermion model |^ . The mirror fermion model is 
discussed in ref. ||]. In this contribution we use 
the reduced staggered formalism QJ^. The basic 
idea here is to couple the two reduced staggered 
flavors to the Higgs field. For this purpose it is 
convenient to introduce the 4x4 matrix fields 

b 

*x = iE(^'^')H(i + £=.+f)x.+6. (1) 



with 7^ • • •74^ ^ (_i)xi+x2+x3+:e4 and 

the sum running over the corners of a hypercube, 
6^ = 0, 1. The one-component staggered fermion 
field X is a real Grassmann variable. In contrast 
to usual staggered fermions (s. eq. (1) of ref. 
we inserted here the factors ^{I—Ex) and ^{l+£x) 
which are used in the reduced staggered formal- 
ism to project the usual staggered fields x ^^J^d x 
to the odd and even sites of the hypercubic lat- 
tice. The form (|l|) implies the following structure 
for the matrix fields 



Vfl 



i^R 





(2) 



where ipi^, ipRj ''pL and iI^r are 2x2 matrices. 
The row (column) indices of the and ^r fields 
act as Weyl-spinor (flavor) labels, and vice versa 
for ■0^ and ipj^. When using the relation (||) and 
introducing the 4x4 matrix for the 0(4) Higgs 
field (A, 



$ = 




(^t 



* Presented by W. Bock 



= -X1'^m7m, (3) 

one can show that the action 

Sf = -^[^iTr(W,7^*,+^-¥,+/.7^*,) 

X /J, 

+ yTr(*,vI;^$^)] (4) 

reduces in the classical continuum limit to the ac- 
tion of the fermion-Higgs sector of the StM with 
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Figure 1. Upper figure: Phase diagram at A = cx) 
with Nd — 2. Lower figure: rria/vR as a function 
of mp/vR. 



one mass-degenerate isospin doublet. After in- 
serting (|l]) and (^) into eq. (||) the final form of 
the fermionic action in terms of the x fields reads 

Sf^-^YI XxXx+ftiVt^x + y£xCt^x'(Pf,a,) 1 (5) 

x^ 

where Ip^^ — J^b Vti,x-b is the average of the 
scalar field over a lattice hypercube and fj^^x — 
(^_lYi+-+x^-i ^ (-^^ ^ (^_lY^+i+-+x^ are the 

usual staggered sign factors. The total form of 
the action is given hy S — Sf + Sh, where Sh — 

SxP't 'Pax^a,x+p, — faxVax — HVaxVax--l)'^] 

is the pure scalar field action. The action S is 
invariant under the so-called staggered fermion 
(SF) symmetry group which includes shifts by one 



lattice distance, 90° rotations, lattice parity and 
the global U(l) symmetry, Xx e^°"^''Xx- This 
invariance of 5* ensures the staggered flavor inter- 
pretation in the scaling region. 

S is, however, not invariant under the full 0(4) 
flavor group: There are two operators with di- 
mension four which are generated by the quan- 
tum fluctuations and which are invariant under 
the SF symmetry group, but break 0(4), 

x^ x^ 

In order to recover the full 0(4) symmetry one 
has in principle to add these operators as coun- 
terterms to the action S S + eoO^^^ + SqO'^^^ 
and tune the coefficients Eq and Sq as a function 
of the bare parameters such that the 0(4) invari- 
ance gets restored in the scaling region. Here we 
shall not add these counterterms to the action. 
However, we will show in the next section that 
the effect of the symmetry breaking is small in 
the parameter region of interest. 

Since we are interested in the largest possible 
renormalized couplings we have fixed in the nu- 
merical simulation A = c». For the use of the Hy- 
brid Monte Carlo algorithm it is necessary to use 
two mass-degenerate doublets, Njj = 2. The K-y 
phase diagram is shown in fig. 1. There are four 
different phases, a paramagnetic (PM), a bro- 
ken or ferromagnetic (FM), an antiferromagnetic 
(AM) and a ferrimagnetic (FI) phase. The var- 
ious symbols mark the points in the FM phase 
where we carried out numerical simulations on 
lattices ranging in size from 6'^24 to 16'^24. 

2. 0(4) SYMMETRY BREAKING 

To estimate the amount of 0(4) symmetry 
breaking we have computed the one- and two- 
point functions in the FM phase using renormal- 
ized perturbation theory. We can decompose the 
scalar field in the FM phase in a Higgs mode, 
tj, and three Goldstone modes, 7r°, a = 1,2,3, 
according to (pR^ = {vr + crR)e'^^ + "^r^Iii where 
{e"} form an orthogonal set of 0(4) unit vec- 
tors, which is arbitrary when neglecting fermion 
loop effects. This arbitrariness is removed after 
taking into account the one fermion loop contri- 
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Figure 2. Goldstone mass as a function of y. 



bution to the vacuum expectation value. One 
can show that the direction of spontaneous 
symmetry breaking is compatible with the 
one fermion loop correction only if the = 

(±1,0,0,0), (0,0,0,±1), (±1,±1,0,0)/V2, 
(±1, ±1, ±1, ±l)/2. From a computation 
of the one fermion loop effective potential for 
small hr one finds that only the directions — 
(±l,±l,±l,±l)/2 correspond to local minima, 
the others are saddle points or local maxima. 

The calculation of the one fermion loop contri- 
bution to the self-energy Yj^^{p) shows that both 
terms in eq. (|^) are generated by quantum fluc- 
tuations and should therefore be included in the 
tree level effective action with JqjSo ^R^eR. 
Explicit expressions for E^^(p) and the coefh- 
cients £r and 5r are derived in rcf. Q in renor- 
malized perturbation theory, er — f^{m.p)NDy'^, 
Sr = fs{'mp)NDyR- The lattice integrals feimp) 
and fsimp), whose explicit form is given in 
ref. 1^, can be computed numerically, fe^mp) = 
0.0054,0.0043, fsimp) = 0.026,0.017 for mp ^ 
0, 0.5. From the one loop result for the renor- 
malized propagator we can read off the following 
estimate for the Goldstone mass. 



1 + Sr/A- 



(7) 



shows that the numerical values (squares) for 
ttItt are very small in the scaling region with 
y w 3.6 — 4.0. Moreover the numerical results 
for TOtt are in good agreement with the analytic 
prediction ^ (diamonds) after inserting the mea- 
sured values for yR and mp. This motivates us 
to take also the corrections for the renormalized 
field expectation value and Higgs mass in eq. (0) 
seriously and to define corrected couplings, y'j^ — 



VRil + ¥)'/^ and A', 



2v; 



-(l-^)(l + ^) = 



The Efl-tcrm in the effective action gives rise to 
the non-zero value of the Goldstone mass. Fig. 2 



A measure for the 0(4) symmetry breaking cor- 
rections is given by the ratios Ry = (yR — y^) /y^ 
and Rx = W2\r - a/2A^)/a/2A'^. A numerical 
calculation of these ratios gives \Ry\ < 5% and 
|i?A| < 7%, in a parameter region with mp < 0.5 
and ma < 0.7, which shows that the symmetry 
breaking effects are small. 



3. RESULTS OF THE SIMULATION 

Since the effect of the symmetry breaking is 
small we have computed the renormalized cou- 
plings from the usual tree level relations yR = 
mp/vR and Xr = m^/2w^, where the renormal- 
ized field expectation value is defined as vr = 
v/y/Z^. Here v is the unrenormalized scalar 
field expectation value and the wave-function 
renormalization constant of the Goldstone prop- 
agator. For the determination of the quantities 
mp, ma- and we have measured the fermion, 
a particle and Goldstone propagators in momen- 
tum space. The fermion propagator could be 
well described for all momenta by a one pole 
Ansatz, which is characteristic for weakly inter- 
acting fermions. For the Goldstone propagator 
we have displayed a typical example in fig. 3. 
The inverse propagator G^^ is plotted here as 
a function of p^, the square of the lattice mo- 
mentum. The numerical data (crosses) exhibit 
a significant curvature at small p'^ which can be 
described by the one fermion loop contribution 
to the self-energy. This non-linear p^ dependence 
can be parametrized by the Ansatz G~\{p) = 
{p"^ +ml^^ + T,sub{p)) /Za,7r, where Ssub is the sub- 
tracted one fermion loop self-energy. The circles 
in fig. 3 were obtained by fitting this Ansatz to 
the numerical data. The fact that the one loop 
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Figure 3. ^(j?) as a function of p^. 



Ansatz is sufficient to describe tfie numerical re- 
sults perfectly over a large momentum interval 
indicates already that the renormalized couplings 
are small. This fitting method allows us to de- 
termine rria and accurately, also on small vol- 
umes. 

As a next step we have to extrapolate the fi- 
nite volume results for and Xr to the infinite 
volume. We carried out simulations on lattices of 
size L^24 with L ranging from 6 to 16. If the spec- 
trum contains massless Goldstone bosons, this 
gives rise to a l/L^ dependence of the finite vol- 
ume quantities. Since the Goldstone particles are 
massive in our model we expect deviations from 
the linear l/Lp' dependence when the volume in- 
creases beyond the Goldstone correlation length, 
L > 0(l/m7r)- The fact that we did not observe 
significant deviations gives further evidence that 
the symmetry breaking effects are small. 

In the lower graph of fig. 1 we display the infi- 
nite volume results for the ratios m„/vR = \/2Xr 
and mpjvB = Vr- The symbols in the upper 
and lower diagrams of fig. 1 match, so that one 
can see where in the phase diagram the results 
for the ratios have been obtained. It can be seen 
that the numerical values for neither ratio change 
when lowering k beyond k = 0, while keeping the 
cut-off roughly constant. The vr values of these 



points vary from 0.08 to 0.27. The arrows in fig. 1 
mark the tree level unitarity bounds for Xr and 
Hr. The graph shows that the points obtained in 
the regions (II) and (III) of the phase diagram 
(see fig. 1) are still very close to these values, 
which indicates that the renormalized couplings 
are not very strong. The solid line encloses the 
allowed regions obtained by integrating the one 
loop /3 functions from infinite couplings at the 
cut-off downward to the renormalization scale. 
The cut-off was adjusted such that the agreement 
with the numerical data is best. It is remarkable 
that the shape is in reasonable agreement with 
our data. Fig. 1 shows that the Yukawa inter- 
action gives a slight increase in Xr. From fig. 1 
we can read off an upper bound for ma/vR and 
fnp/vR: For m^'^Q.l /a, we find rricr/vR^A and 
mp /vr^2.6. From experience in the 0(4) model 
with various regularizations, we expect that these 
numbers for the upper bounds may be stretched 
by perhaps 20-30%. 

All in all we conclude that the renormalized 
quartic and Yukawa couplings are in accordance 
with triviality and that they cannot be strong, 
unless the cut-off is unacceptably low. 

The numerical calculations were performed on 
the CRAY Y-MP4/464 at SARA, Amsterdam, on 
the S600 at RWTH Aachen and on the CRAY Y- 
MP/832 at HLRZ Juhch. This research was sup- 
ported by the "Stichting voor Fundamentcel On- 
derzoek der Materie (FOM)" and by the "Sticht- 
ing Nationale Computer Faciliteiten (NCF)" . 
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